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ABSTRACT 

Compact remnants - stellar mass black holes and neutron stars formed in the inner few 
parsec of galactic centres are predicted to sink into the central pa rsec due to dyn amical 
friction on low mass stars, forming a high concentration cusp (jMorris. Tool . Same 
physical region may also contain very high density molecular clouds and accretion 
discs that are needed to fuel SMBH activity. Here we estimate gas capture rates 
onto the cusp of stellar remnants, and the resulting X-ray luminosity, as a function 
of the accretion disc mass. At low disc masses, most compact objects are too dim 
to be observable, whereas in the high disc case most of them are accreting at their 
Eddington rates. We find that for low accretion disc masses, compact remnant cusps 
may be more luminous than the central SMBHs. This "diffuse" emission may be of 
importance for local moderately bright AGN, especially Low Luminosity AGN. We 
also briefly discuss how this expected emission can be used to put constraints on the 
black hole cusp near our Galactic Centre. 
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1 INTRODUCTION 

Super-massive black holes are predicted to be sur- 
rounded by clusters of compact stellar remnants, i.e., 
stellar mass black holes and n eutron stars. Dynami- 
cal friction (jChandrasekhar. 19431 ) on less massive back- 
ground stars transfers these remnants closer to the SMBH 
l|Morris. 19911 ). For the case of our Galactic Centre, stel- 
lar mass black holes can be collected into the cen- 
tral parsec from a sphere with radius about 4 par- 
sec l|Morris. 19931 ). iMiralda-Escude fc Gould (2000l ) up- 
dated and quantified Morris (19931 ) estimates and showed 
that black holes will force most of low mass stars out 
of the central region with size f?bh ~ 0.7 parsec around 
Sgr A*. In their model, black holes form a central clus- 
ter of this size and number about 20, 000. The total 
mass of these black holes is around 5 % of Sgr A*'s mass 
ISchodel et al.. 200jGhez et al.. 20031 ). iFreitag et al. (20061 ) 
have recently modelled the dynamical mass segregation in 
the central few parsec of our Galaxy with state of the 
art Monte C arlo methods, largely confirming re sults of 
iMorris (199t ) and iMiralda-Escude fc Gould (2000D . In ad- 
dition, it was shown that the central parsec contains neutron 
stars in numbers comparable to that of the black holes. 
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Recent observations of young massive stars near 
Sgr A* suggest that these stars may also be born 

in situ with an unusually top-heavy mass function 

(|Paumard et al.. 200dNavakshin fc Sunvaev. 200jNavakshin. 2006 
The remnants of the young stars born in situ will contribute 
to the cusp of black holes and neutron stars formed by the 
dynamical friction. 

iMorris (19931 ) discussed the possibility that compact 
remnants will be accreting gas from molecular clouds that 
may be unusually dense near our GC. In this paper we set 
out to explore this idea in more detail and also speculate 
on the implications for observations of AGN in general. Our 
model includes a realistic form for the radiative efficiency of 
these sources motivated by the observations of radiatively 
inefficient accretion flows in X-ray binaries. We integrate 
over the black hole velocity distribution to obtain a time- 
average estimate of total X-ray luminosity of the cluster as 
a function of the mass of the embedded gaseous accretion 
disc. 

We find that at low disc masses, most of the stellar 
remnants accrete at highly sub-Eddington radiatively ineffi- 
cient rates, and are quite dim. However, there are always few 
that almost co-move with the gas in the disc. These objects 
achieve relatively large accretion rates. Since the SMBH it- 
self becomes radiatively inefficient at low disc masses, our 
model predicts that clusters of stellar remnants may ac- 
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tually outshine the SMBH in this Low Luminosity AGN 
(LLAGN) regime. At the high luminosity, high disc mass 
regime, the super-massive black hole should dominate bolo- 
metrically. Stellar remnants may however be important in 
certain wavelengths. 



2 THE MODEL 

2.1 Bondi-Hoyle gas capture rates 

We shall first make order of magnitude estimates to motivate 
a more detailed treatment. A stellar mass black hole of mass 
Mbh, ploughing throu gh a gas cloud, is capturing gas at the 
iBondi fc Hovle fl944h accretion rate: 
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where p, Cs and Av are the ambient gas density, sound speed 
and the relative velocity between the black hole and the gas, 
respectively. 

Consider first the situation in the central parsec of our 
Galaxy. Stellar velocity dispersion, a, is of the order of a few 
hundred km/sec there. The number of black holes, AA^'bh, 
moving with respect to gas with a relative velocity smaller 
than All is of the order of ATVbh ~ Nhh{A.v /veec)'' fv , where 
A^'bh 2> 1 is the total number of stellar mass black holes, 
Vcsc ~ (T is the escape velocity, and <C 1 is the vol- 
ume filling factor of cold gas in the region. For an estimate, 
take fv = 0.01 and A^bh = lO". In that case, AA^bh ~ 1 
for Av ^ 50 km/sec. There thus should be few stellar mass 
black holes travelling through the gas at At; ~ 50 km/sec. 
Neglecting the gas sound speed in equation 1, the Bondi- 
Hoyle capture rate is 

Mcapt ~ 4.8 X lO-'"* Moyear"' 



(2) 



where no is the gas density in Hydrogen atoms per cm'', 
p is mean molecular weight in units of proton mass, and 
Av5o = A«/(50) km/sec. 

This accretion rate is quite small compared to the Ed- 
dington limit of MEdd « 2 x 10"'^( Mbh/10 Mq) year"^ 
if no = 1. However, gas density can be much higher. The 
hot X-ray emitting gas in the inner parsec of our Galax y 
has density of n ~ 10 - 100 cm'^ jBaganofr et al.. 2003h . 
If we assume cool gas clouds to be in pressure equilib- 
rium with the hot gas, then the cool gas density could be 
as high as 10^ — 10® cm"''. Therefore we could in princi- 
ple expect several stellar mass black holes accreting gas at 



rates corresponding to 
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10"^ of their Eddington 
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where pcrit is the critical density at which the disk is 
marginally self-gravitating. Numerically, 

Ms 



n = 



— = 2.8 X 10^^cm"^(: 
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where Ms = 10-*Msmbh/ Mq, iJo.i = ii/0.1 parsec and 
rUp is proton's mass. For such large gas densities, even 
black holes with very large differential velocity Av, e.g., 
Av ^ 1000 km/s, would accrete at their super-Eddington 
rates: 
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In practice, this implies that every black hole crossing the 
AGN accretion disc captures the material at or above its 
Eddington accretion rate. 

2.2 Maximum accretion rate 

A black hole co-moving with the gas, e.g. with At; = 0, is 
predicted to accrete at a very large accretion rate by the 
equation 1. However, in reality gas must be in a Keplerian 
differential rotation about the SMBH. Thus, even if gas ex- 
actly coincident with the stellar mass black hole has zero 
relative velocity, the neighbouring regions will have Av 7^ 0. 
In such situation, the sphere of infiuence of the stellar mass 
black hole is th e Hill's radius rn = -R( Mbh/SMsMBu)^''^ 
l|Lissauer. 1987h . and the maximum accretion rate is of or- 
der of 



Mh = ^TTThPCs ■ 

We thus require 
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rates. Since Sgr A* pr oduces only ~ 10 erg/sec in X-rays 
(|Baganoff et al.. 20031 ). the brighter of the accreting black 
holes can actually outshine Sgr A* in X-rays. 

In the case of a reasonably bright quasar or 
an AGN accreting through a parsec scale accretion 
disc, we can estimate the gas density in the disc 
in the following way. In order to be stable against 
self-gravitational collapse, the disc's iToomre (19641 ) 
Q-parameter should be smaller than unity (e.g., 
iGammie. 200llGoodman. 200^hompson et al.. 200dNavakshin. 
The gas density in the disc satisfies 

-1 MsMBH 



2.3 Expected accretion luminosity 

As is well known, standard accr etion discs have a very hig h 
radiative efficiency: eo = O.lc^ (jShakura fc Sunyaev. 19731 ). 
However, at low dimensionless accretion rates, rh = 
M / MEdd ^ nicrit, radiative efficiency is predicted to be 
much sm aller since electrons i n the flow may be cooler than 
protons (|Naravan fc Yi. 19951) and since ga s outflows may 
be important (iBlandfo rd fc Begelman. _ 1999l) . For definitive- 
ness, we follow the iNaravan fc Yi (19951 ) model, in which 
e — eo(ni/rhcrit) for rh < rhcrit, and rhcrit = 0.01. This 
mode l seems to be c ompatible with X-ray observations 
(e.g., lEsin et al.. 199?! ). Joining the two accretion regimes 
smoothly, and assuming that roughly 10% of the bolometric 
luminosity comes out in the X-ray band, we write 



e(m) 



m 



(8) 



m -I- merit 

where eo,x = 0.01. For accreting neutron stars, the radiative 
efficiency is assumed to be always large and equal to eo,x 
due to the presence of hard surface. 



Accretion and emission from the black hole 
cluster 



Black holes in coUisionless cusps around SMBHs are ex- 
pected to have the isotropic iBahcall fc Wolf (19761 ) space- 
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velocity distribution. This distribution is a good approx- 
imation to the rec ent detailed Monte-Carlo simulations 
ijFreitag et al.. 20061 ). in which the radial black hole den- 
sity distribution found to follow the law pbh{R) oc f?~^'^*. 
To keep our analysis analytically simple, we approximate 
this law to pbh{R) oc R~^^^, in which case the black hole 
space and velocity distribution is particularly simple (e.g., 
iMiralda-Escude fc Gould. 2000l ): 



•' ^ ' ' 1 0, otherwise. 



(9) 



where we also assume R < Rhh, the radius of the black 
hole cluster around the SMBH. The constant Co is set by 
requiring A''bh ~ J dH J civ/(R, v). Using this velocity dis- 
tribution, we calculate the total accretion rate onto the stel- 
lar mass black holes, assuming that gas moves on circular 
Keplerian orbits in a thin disk of half thickness H <^ R. 



Mtot = / dRp(R) / dv/(R,v) 



'p(R) 



(10) 



Note that, except for very high accretion rates, Mbh = 
Afcapt oc p, and hence the expression under the velocity in- 
tegral does not depend on p(R). We make one further ap- 
proximation, writing VeBc{R) ~ fcsc(i?bh)- This is justified 
since most emission arises from the edge of the black hole 
cluster, i.e. R ~ -Rbh- This is so because that region contains 
most of black holes and they also have the smallest Av, as 
Vesc is smaller than in the inner regions. We thus have. 



Mu 



dvCo 



c(flbh) 



Mbh 
p{Rhh, 



(11) 



In this simple model, the total accretion rate onto stellar 
mass black holes scales linearly with the disc mass. The total 
luminosity of the black hole cluster is given by a similar 
equation to [10] but with the additional e(m) (equation [8]) 
included under the sign of the integral. 

A more detailed treatment of this problem is under- 
way (Deegan & Nayakshin, in preparation). In this Monte 
Carlo-like approach, we follow trajectories of accreting ob- 
jects explicitly, and treat both capture of gas and "small" 
scale disc accretion in a time-dependent manner. The results 
are consistent with those obtained in this paper but show a 
large time-variability, as expected. 

2.5 Model of disk accretion for the SMBH 

In the absence of a reliable theory which would connect ac- 
cretion disk properties at parsec scales to the inner regions of 
accretion flows on SM BH at both low and high mass accre - 
tion rate regimes fe.g.. lGoodman. 200jCuadra et al.. 20051 ). 
we model the accretion onto the SMBH in the follow- 
ing way. We estimate the viscous time scale for ac- 
cretion through the disk as tviac = {H/ R)~'^Q~^ 
IShakura &: Sunvaev. 197a ). The accretion rate onto the 
SMBH is then 
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H 
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The disc height scale H is calculated assuming the hydro- 
static balance as long as the disc is non-self-gravitating 
(i.e., relatively low-mass). For simplicity we assume 
disk temperature is equal to 100 Kelvin. The results 



depend weakly on this value as H/R oc in this 

regime. For more massive self-gravitating disks (see, e.g., 

iPaczvriski. 197dGammie. 200llGoodman. 200dNavakshin fc Cuadra. 20051 ). 
we assume Q ~ 1, which translates into H/R ~ 
A/disc/MsMBH. The SMBH luminosity is calculated in 
the same way as for the stellar mass black holes, i.e. 
isMBH = e{m)MsMBH, where e(m) is given by equation |8l 



3 RESULTS AND OBSERVATIONAL 
IMPLICATIONS 

With this model we now calculate two examples. The first 
one is based on th e up to date models of the central par- 
sec of our Galaxy IjFreitag et al.. 20061 . e.g.,), and the sec- 
ond is for the M = 10* M© SMBH case. For both of these 
we assumed the total mass of the stellar mass black holes 
to be 5% of the mass of the central SMBH. For the GC, 
this yields iVbh = 24,000, whereas iVbh = 720,000 for 
the more massive galactic nucleus. The black hole cluster 
is assum ed to extend to Rhh = 0.7 pc for Sgr A* case 
( Miralda -Escude fc Gould. 2000l ). and to i?bh = 2 pc for the 
other case. We also calculate the expected X-ray emission 
from neutron stars in the clus ter, taken to be as nu merous as 
the stellar mass black holes (|Freitag et al.. 20061 ) but have 
the mass of 1.4 Mq. Note that in this paper we consider 
neutron stars to be magnetic field free. If their fields are 
strong enough, the "propell er effect" may reduc e the X-ray 
luminosity of neutron stars (|Menou et al.. 19991 ). 

Figure [T] shows the X-ray luminosities of: the SMBH 
(dotted), the cluster of stellar mass black holes (solid), and 
the neutron star cluster (dashed) as a function of the accre- 
tion disc mass. The left panel shows the GC case, whereas 
the right panel shows the Msmbh = 10** M© SMBH case. 
The main prediction of Figure [T] is that at high disc masses 
the black hole cluster emission is negligible compared to that 
of the SMBH, but at smaller disc masses the situation is re- 
versed. This is understood as following. At the highest disc 
masses, all of the components considered accrete at rates ap- 
proaching their respective Eddington limits. As the SMBH 
is more massive than the cusps we model, it is brighter than 
the stellar mass compact remnants. At very low disc masses, 
on the contrary, the SMBH accretes in the radiatively inef- 
ficient regime, and it is quite dim. At the same time, the 
accretion rates onto the stellar mass black holes span a wide 
range according to a range in the relative velocities Ad. 
Though small in numbers, black holes with small relative 
velocities dominate the radiative output of the cluster since 
the radiative efficiency for them is largest (see equation [S} . 
We now discuss the observational implications of our results 
in greater detail. 



3.1 "Fake X-ray binaries" in Sgr A* 

It is only in the centre of our own Galaxy that we have ob- 
servational constraints on the distribution of gas and X-ray 
sources in the inner parsec. The dominant reservoir of cold 
gas i n the inner parsec of our GC is the so-c alled minis- 
prial (IScoviUe et al.. 200dPaumard et al.. 20041 ). This mor- 
phologically complicated feature covers a good fraction of 
the inner parsec in projection. It is believed to be mildly 
geometrically thick, and have temperature and density of 
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10' 10^ 10^ lO'* 10^ 10' 10"^ 10' 10*" 

Disk mass Disk mass 

Figure 1. The X-ray luminosity (in erg s^^) of the different accreting components versus the accretion disc mass (in units of Mq). 
In particular, the solid, dashed and dotted curves show the luminosities of the black hole cusp, the neutron star cusp, and the SMBH, 
respectively. The left panel is for the case of Sgr A*, and the right one for a SMBH with Mj^jj = 10* Mq. Note that at low disc masses 
(low SMBH luminosities), the cusps of compact remnants are brighter than the central massive black hole. 



order of T ~ 6 x 10^ K and no ^ lO*, respectively. The high 
temperature is maintained by photo-ionisation via radiation 
of nearby luminous massive star s. The mass of the m ini- 
spiral was estimated at ~ 50 M© (|Paumard et al.. 20041 ). al- 
though observational uncertainties are large, and the real 
mass could be as high a factor of 10 larger (T. Paumard, 
private communication) . 

With Mdisc = 50 M© , the left panel of Figure 1 predicts 
that the X-ray luminosity of Sgr A* should be ^fewxlO"^* 
erg/s. This is about a factor of ten larger than actually ob- 
served (|Baganoff et al.. 20031 ). Given many approximations 
and uncertainties that are used in our simple model, the dis- 
agreement is not surprising, especially since the fraction of 
radiation power coming out in X-rays is exceptionally low 
for Sgr A* (which is not taken into account in our model). 



3.2 Quasars 

The right panel of Figure 1 shows that the emission from 
the cusp of stellar mass remnants is bolometrically unim- 
portant for a quasar accreting at a good fraction of its Ed- 
dington limit. As already discussed, the result is not surpris- 
ing, as the total mass of the stellar mass black hole cluster 
is ~ 5% of the SMBH mass in this model. However, the 
spectrum of stellar mass black holes may be sufficiently dif- 
ferent from that of the quasar that the former may still 
contribute to some wavelength range. In particular, the ac- 
cretion disks of stellar mass black holes are systematically 
hotter than those of much more massive quasar black holes 
(jShakura fc Sunyaev. 1973 ). The black hole cusp could then 
contribute to the soft X-ray band emission of quasars where 
quasar disks are not expected to radiate much energy (ex- 
cept from a corona or a jet). 



Figure 1 further shows that the population of 24,000 
black holes would emit around 10"^^ erg/s in X-rays. Most of 
this luminosity comes from a few to a dozen of X-ray sources 
with X-ray luminosities ranging between 10^^ to 10^^ erg/s. 
It thus seems possible that some of the unusual X-ray sources 
in the inner parsec detected by iMuno et al. (20051 ) are the 
isolated accreting black holes or neutron stars masquerading 
themselves as X-ray binaries. Such compact isolated accre- 
tors will probably have unusual time-variability properties 
as their discs may be much larger than the typical discs of 
X-ray binaries, and indeed they are missing the perturbing 
influence of the secondary. On the other hand, accretion disc 
feeding in these sources will be variable itself, leading to vari- 
ability on variety of time scales. Illumination of nearby gas 
clouds by X-rays may also lead to reprocessed infra-red, op- 
tical and UV emission. Finally, the sources may leave trails 
of denser and likely hotter gas behind them as they plough 
through the gas. 



3.3 LLAGN 

Low Luminosity AGN are the systems where the accreting 
stellar mass black holes may potentially complicate inter- 
pretation of the data, as the X-ray luminosities of these sys- 
tems roughly range from Lx ~ 10''^ erg/s to Lx ~ 10*^. 
Further, unlike Sgr A*, emission of these black holes is spa- 
tially unresolvable from that of the LLAGN. It is desirable 
to the broad band spectrum of the black hole cusp in the 
future work. This may be non trivial if some of the bright 
systems are in poorly understood super-Eddington accre- 
tion rate re gime, similar to the famous X-ray binary system 
SS443 (e.g.. lFabrika. 200d ). 

A way to distinguish between the SMBH and black hole 
cusp emission is via time-dependent observations. The emis- 
sion from a relativistic region of a SMBH might vary on 
time scales of hours to days, as observed in both Sgr A* 
llBaganoff et al.. 200ll ). and in more distant AGN (e.g., 
lEdelson et al.. 20021 ). The emission of a large number of sta- 
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tistically independent black holes and neutron stars should Thompson T. A., Quataert E., Murray N., 2005, ApJ, 630, 167 

be considerably less variable. Toomre A., 1964, ApJ, 139, 1217 



4 CONCLUSIONS 

The presence of stellar mass black holes 
and neutron stars in a cusp around a 
SMBH is a robust theoretical prediction (e.g., 
iMorris. lOOjMiralda-Escude fc Gould. 200dFreitag et al.. 200^ 1. 
The underlying physics - dynamical friction of these stellar 
remnants on a background stellar population - is very well 
understood. Here we attempted to estimate the accretion 
rates with which these objects would be capturing gas 
from a massive accretion disc. We also estimated the 
X-ray luminosities of these objects. We found that the 
total luminosity of stellar mass remnant cusps might be 
comparable to or larger than the luminosity of the central 
SMBH at low gas densities. These conditions are plausibly 
met in LLAGN and the centre of our Galaxy. Observational 
constraints on emission from these sources could be used to 
constrain population of these sources in the central parsecs 
of galaxies. 
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